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The results of measurements of "fN total cross section, obtained by 
the method of photoproduction processes registration at the Baksan 
Underground Scintillation Telescope, are presented. These data at 
energies ^fs = 40 — 130 GeV confirm the effect of more rapid photon- 
hadron cross-section rise as compared to the hadron-hadron ones. It 
is shown, the increasing of the additive quark number in the products 
of photon hadronization can be one of the causes responsible for this 
effect. On the basis of the analysis of experimental data on both 7^ 
and 77 total cross sections, the status of direct and indirect cross- 
section measurements is discussed. 
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1 Introduction. Status of the data. 



Since the first measurements from ZEUS |T| using 1995 data, information 
about the proton (nucleon) structure function F 2 at low Q 2 and cross 
section <7 7 tv(s) at high energies generates a lot of interest. In the energy range 
y/s = 100-^250 GeV, the total cross section of 7p interaction were determined 
by indirect method || on the basis of DESY data on proton structure func- 
tions |IJ. At several points of energy scale around \fs = 200 GeV, the cross 
section was measured by direct method — by photoproduction processes reg- 
istration H |J. The results obtained testify to more rapid photon-hadron 
cross section rise in comparison with the hadron-hadron ones. The physical 
nature of this more rapid rise is the main problem of 'yN physics, which 
is discussed in many papers. Modern interpretation in the framework of a 
Regge parametrization is done in |J. The reviews of experimental data and 
theoretical approaches are given in || [/J. 

In this paper we present the results of measurements of ct 7 at(s) in the 
region yfs = 40^-130 GeV obtained at the Baksan Underground Scintillation 
Telescope (BUST) of Baksan Neutrino Observatory. In this experiment the 
hadronic and electromagnetic cascades produced by cosmic ray muons are 
registered. The BUST is situated at an effective depth of 8.5 x 10 4 g/cm 2 
underground || . So only muons (with energy > 220 GeV) and neutrinos 
can reach the telescope. That is the experiment conditions correspond to 
the pure muon beam with the known energy spectrum. Preliminary results 
were presented earlier ||, however the data on cross section o" 7 at(s) were not 
reported in ||. The data obtained in this paper confirm the effect of more 
rapid rise of photon-hadron cross section. 

The theoretical basis of how to find the 'yN cross section from the data 



on muon interaction with nuclei was developed in [fL0|l , |TTJ. It was shown 
in this Refs., that in BUST experiment the number of hadronic cascades is 
proportional to a^ N (s). 

To interpret the experimental data, we propose a phenomenological the- 
ory based on the Vector Dominance Model (VDM) and the Additive Quark 
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Model (AQM). In this theory the violation of photon-hadron scaling is 
parametrised by two effects: 1) the rise of hadronization probability; 2) the 
rise of hadronization product, which is due to vector resonances decay on 
pions. It is shown, the experimental data on 7 AT cross section allows to find 
the measures of these effects and than to predict the 77 cross section with 
high enough accuracy. 

Let us introduce a quantitative characteristic of the effect of the photon- 
hadron scaling violation. According to the VDM 

(T w (s) =Y / P i^v(s)cr vp (s), (1) 

v 

where P 7 _^y( s ) is the probability of photon conversion to a vector meson V; 
a vp (s) is the total cross section of Vp interaction. Making use the AQM [T2| 
and U(3) symmetry of meson-nucleon interactions gives 

2 

<7v P (s) = t^np (3s/2), V = p, u, 0. (2) 
This allows to separate in ([!]) the contributions of p, u, <fi mesons: 



a 



$(s) = P 7 ^(s)-a Np (3s/2) (3) 



where 



P^^puxpi^) — ^ ] Py— fjvp — ^S^PP ®PP ^np &np)- (4) 

V=p,u>,cf> 

The argument in the right-handed part of Eq. is multiplicatively trans- 
formed (3s/2 instead of s), because the energies of interacting quark pairs 
in Vp and Np systems are different. The photon hadronization probability 
Py-*puj<p ~ P^pujrp = 1/250 is known from the data on 7 —>■ V transitions 
on the mass shells of vector mesons. So p,co,4> contributions (i.e. a^°J(s) 
function) can be expressed only in terms of the data on cross sections of 
nucleon-proton interactions. The function cx^(s), by virtue the specifying it 
experimental data are independent and the physical sense is clear, is referred 
to as the calibration curve for the cross section of 7p interaction. It is clear, 
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the deviation of measured cr 7p (s) from the calibration curve contains an ad- 
ditional information about photon hadronization, which is not contained in 

For constructing the calibration curve in the energy region yfs = 5.93 -r- 
22.97 GeV, we use data from [0, 0. In the region = 30.4 + 62.7 GeV, 
the data of [|l^, [U| on u pp (s) and cr pp (s) are used (it seems to be reasonable to 
assume, that at such energies a pp = <r np , a pp = a pn with a good accuracy); for 
energies ^/s = 200, 546, 900, 1800 GeV we use the pp data from |TJ, [18], 0, |20| . 
These data, subjected to the scale transformation 3s Np /2 = s lp = s, are fitted 
by the calibration curve 

a%(s) = C\n^ + A(l + -^ r ) (5) 

So S + A 02 

C = 0.5777 /ib, ^= 2.198 GeV, A = 95.65 /xb, 
Aoi = 2.774 GeV, A 02 = 3.589 GeV. 

The data and the calibration curve mentioned above, and also various 
experimental points on u 1P (s) (including those of BUST) and two variants of 
their fit, corresponding to the two variants of the high energy points choice, 
are presented in Fig.l. The upper solid line represents the cross section 
a 7P (s) as if the photon hadronizates only in the strongly bound 7r + 7r - system; 
physical bases for invoking of this curve are discussed in section 3. 

One can see from Fig.l, at low energies the photon-proton cross section 
asymptotically approaches the calibration curve. This testifies, that con- 
siderations used at constructing the calibration curve adequately reflect the 
physics at these energies. Essential deviations from the calibration curve 
emerge at ^Js > 8 GeV. At */s = 200 GeV the deviation of the cross sec- 
tion obtained by the photoproduction (direct) method achieves 25%. That 
derived from the extrapolation of ZEUS data |j] (indirect) is even larger. 
The BUST data, obtained by the photoproduction method, are combined and 
fitted with both direct and indirect measurements. The results are shown in 
Fig.l as the dashed and dotted lines correspondingly. 
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* experimental data and 
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Figure 1: Calibration curves and cross section of 7p interaction. 

2 Experimental technique and the results at 
the BUST. 

2.1 The measured variables. 

In 1983 - 1989 at the BUST, the experiment on studying cross section 
da^lE^, E c ,q 2 )/dE c of inelastic scattering of Cosmic Ray (CR) muons on 
atomic nuclei was carried out |9|, [H], 22]. Here E c = — E' is the energy 



transferred from muon into the hadronic cascade, E,. and E' are muon ener- 
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gies in laboratory system before and after interaction. The average number 
of nucleons in the nuclei of the target was estimated as A = 26. The events 
characterized by E c > 700 GeV were selected. 

The average value of the squared transferred four-momentum in our ex- 
periment, — q 2 = Q 2 = 0.5 GeV 2 , is much less than the characteristic squared 
energy in the centre-of-mass system of photon and nucleon, s 7iV > (40 GeV) 2 
(Ery > 700 GeV). This allows us to pose the problem on experimental mea- 
surements of interaction cross section cx 7iV (s) of a real photon with a nucleon. 

The problem of the measurement a~ N in the processes with Q 2 = 0.5 GeV 
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is nontrivial. In |T(| [11] it was shown, that the differential cross section of 
inelastic \iA scattering at small Q can be represented as a product of 
and a certain function. This function depends on Q 2 and other parameters 
taking into account the effects of shadowing in nuclei; its explicit form is given 
in the Refs. mentioned above. Notice, the factorization of (Try JY IS Cl general 
property of fiA cross section at small enough Q 2 ; the question is what is the 
particular value of Q 2 that the factorization can be actually performed at. 
In [TI| |TTJ it is shown that this can be done if Q 2 < 10 GeV 2 . This structure 
of fj, A scattering cross section allows at Q 2 <C 10 GeV 2 to determine the a 7JV 
cross section, in spite of the fact that in each separate act of measurement 
the processes with Q 2 ^ are fixed. In the expression for cross sections of 
inelastic \iA scattering (at E c > 700 GeV inelastic fiA scattering initiates 
a hadronic cascade), it is possible to factorize the dependence from Q 2 , to 
integrate over this dependence and then to express the cross section directly 
through a jN (E c ) [0: 

dafl(E^E c ) _ a AE C 



dE c 2vr E 2 

A 4 



a 7N (E c )-F A (E^,E c ), (6) 



where F^E^, E c ) is a function which takes into account shadowing of nucle- 
ons in a nucleus. 

In our experiment Q 2 is determined as 



fo(E c y J " y c/ Je c dE^ JQmin dEJQ 
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where N^E^) is the integral energy spectrum of CR muons (see Eq. (0) )• At 
calculating the integral over Q 2 , we use the expression of d 2 a(E^, E c , Q 2 ) / dE c dQ 2 
and g 1 n(E c ) from (this can be done because a 1 N{E c ) describes the ex- 
perimental data fairly well) and restricted the integration range from Q 2 min 
to 10 GeV 2 . The contribution of the range Q 2 > 10 Gev 2 in Q 2 value is 
< 1%. 

The integral energy spectrum of muons for zenith angles 9 < 25° was 
obtained in the special gauge experiment p2) and with an accuracy of 5% 
approximated by the expression: 

N,(> E ] 9,0) = k-(E + E H )- 2 - 83 , (7) 

where E H = 230 GeV, k = const. At 9 > 25° we used the calculated muon 



spectrum obtained in Ref. p3| 



The BUST has sizes 16 x 16 x 11 m 3 and consists of four horizontal and four 



vertical scintillation planes |J. In experiment |9], [H], [22| the four horizontal 
scintillation planes of the BUST together with the concrete layers between 
them were used as a calorimeter that registered electromagnetic and hadronic 
cascades generated by CR muons either in the rock immediately above the 
telescope (A = 25) or in the matter of the telescope (A = 27). Experimen- 
tally measured values are the differential spectra of energy deposition in the 
calorimeter from the electromagnetic (j = e) and hadronic (j = h) cascades 

JN i( e) _ f 3N,(E r , 0, 4.) d^>(E„, E e )_ ^ ^ ^ ^ ^ 



de J dE^ dE c 

(8) 

Here e is the sum of energy depositions in four horizontal scintillation layers 
of the BUST; 

dN^E/j,; 9, cfy/dE^ is the differential spectrum of muons in a given direction 
at the BUST location; 

da^\{E^, E c ) j dE c is the production cross section of electromagnetic or hadronic 
cascade with energy E c by muon with energy E^ 

n is the number of nuclei in a unit volume of the detector (or of the rock); 
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W(E C , e; 9, 0) is the probability of that a cascade of energy E c deposits the 
energy e in the telescope; 

S(9, <p) is the area of the telescope in a direction 9, (p. 

At calculating the probability W(E C , e; 9, 4>), electromagnetic and hadronic 
cascade curves, obtained by Monte Carlo simulation in view of real structure 
of the facility, are used. W(E c ,e;9,<p) takes into account the integration 
over the thickness of the target, therefore only the integration over S(9, 0) 
remains in Eq. (||). 

Muons produce cascades in the following processes: 

a) bremsstrahlung of photons, 

b) production of e + e~ - pairs, 

c) production of high-energy knock-on electrons, 

d) inelastic scattering on nuclei. 

The hadronic cascades are produced only in the process (d), the processes 
(a), (b), (c) produce the electromagnetic cascades. For the cross sections of 
the processes (a), (b), (c) the expressions were taken from Refs. [23, ^6 . 



To obtain the cross section a^ N (E c ) we use the ratio of numbers of 
hadronic and electromagnetic cascades 

having assumed the cross sections of electromagnetic processes are known. 
Making use of the ratio R(Ae) has an important advantage in comparison 
with the procedure of <t 7 at determination directly from the data on hadronic 
cascades. The possible uncertainties in determination of A^(> E;9,<j)), 
W(E c ,e;9,<p), S(9,<f)) (see Eq.(|8|)) are nearly completely disappear in the 
ratio @. 

In an individual event we do not try to restore the cascade energy E c on 
the energy deposition e, because the error is equal to 30 — 60%. However, 
with a good enough accuracy of 2 — 3%, it is possible to calculate the av- 
erage energies of electromagnetic and hadronic cascades creating the energy 
deposition in the facility in a given interval Si < e < £ i+ %. 
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The properties of the BUST allow us to distinguish cascades generated on 
different muon trajectories (in the cases when the cascade is accompanied by 
a muon group), and two consecutive cascades on one and the same trajectory 
and, thus, to exclude errors in the determination of e caused by these factors. 



2.2 Separation of cascades. 

To separate hadronic and electromagnetic cascades, we used the number of 
7r — > fi — > e decays recorded in the event @, 21]. In hadronic cascades the 
number of 7r — /i — e decays is 20-25 times more, than in electromagnetic 
ones. Recording the delayed pulse on the screen of 10-beams oscillograph is 
the feature of /jl — e decay (each beam corresponds to the plane of the BUST). 
The separation procedure is based on the results of electromagnetic cascades 
modeling and consists in the following. 

1) The probability P(m, e) of recording of m decays in the electromagnetic 
cascade, which creates the energy deposition e in the facility, is calculated 



21 



2) The separation criterion mo(e) is determined, such that the probability 
X of recording the number of decays m > mo(e) in the electromagnetic 
cascade is < 10~ 2 , 



'»() 



J2 P(m;e) = 1 - x(m ,e) - 0.99. (10) 



m=0 



3) Cascades with m < mo(e) are considered as electromagnetic ones. 

4) The calculated number of electromagnetic cascades with m > mo(e) 
is subtracted from the number of cascades with m > m (e) (i.e. it is taken 
into account that the "tails" of distributions P(m,e) fall into the area m > 
m (e)). 

On the basis of (0) one can write 

N(m < m {e)) = N e (e) - N e (e) ■ X {m ,e), 
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where N e (e) is the actual number of electromagnetic cascades. Therefore the 
total number of electromagnetic cascades in the interval Ae is 



The remaining cascades are considered as hadronic ones. 

5) We do not take into account hadronic cascades falling into the area of 
electromagnetic ones. The calculations show, that the fraction of hadronic 
cascades with m < m (e) does not exceed 3 %. It is much less than the 
statistical error, which is estimated as y/Nh- 

Thus, the efficiency of the electromagnetic cascades separation from the 
hadronic ones is > 99%; in this case the efficiency of separation of hadronic 
cascades from electromagnetic ones is ~ 97%. 

2.3 The results of measurements. 

In Table 1 the data are presented on the numbers of electromagnetic N e (Ae) 
and hadronic Nh(Ae) cascades recorded in various intervals of energy depo- 
sition in the facility Ae = e m in + e max . The last interval of Ae is not limited 
from the above: e > 573 GeV. The spectra are nonmonotone, because at 
small e the observation time was smaller (the total time was T rec = 31650 
hours) . 

The average energy E^ of hadronic cascades, creating energy deposi- 
tion in an interval e m i n < e < e max , is identified with the energy of photon, 
E y = E^\ For a given Ae, the photon energy interval E^ min ) 4- -E 7 ( max ) 
is calculated. The last two columns of Table 1 comprise the photon-proton 
energies in the centre-of-mass system and the corresponding total cross sec- 
tions; energies being given in GeV, cross sections being given in fib. 

The information about <t 7 at(s) is extracted from the ratio of numbers of 
hadronic and electromagnetic cascades, as it has mentioned in section 2.1 (Eq. 
(H)). The final results for a lN (E 1 ) are obtained by integration of photon- 
nucleon interaction cross section over the energy interval AE^ = f{Ae) in 
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Table 1: The data on cascades and cross sections of 7iV interaction 
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N e (Ae) 


N h (Ae) 




E-y( ma x) 


E~ 


V* 


ct 7 at 


77 


133 


1050 


96 


744 


1260 


1100 


45.4 


128 ± 14 


133 


230 


1006 


102 


1260 


2100 


1550 


53.9 


146 ± 15 


230 


573 


445 


45 


2100 


4980 


3030 


75.4 


148 ± 23 


573 




36 


4 


4980 




9500 


133.5 


152 ± 80 



view of the cross section dependence on energy. For this purpose the forecast 
of a jN (E y ) offered in |llj and consecutively refining fit of experimental data 
are used, on which the calculated value R ca i{Aei) is determined. The devia- 
tion of experimental value R exp (Aei) from the calculated one determines the 
deviation of actually measured cross section from the forecast. 

Now we can compare the method of a 7 jv measurements used in this ex- 
periment with those used in DESY experiments. 

1. Direct method || |J is based on photoproduction registration. The 
photon being on-mass-shell (the virtuality is negligible) is ensured by strict 
selection of the final hadronic states. 

2. In the indirect method [Q, @] electron in the inclusive process ep — > e+X 
is registered at small Q 2 ; the value of ct^n is obtained by model-dependent 
extrapolation of the structure functions to Q 2 = 0. 

3. In the method described in this paper the photoproduction of slightly 
virtual photon is registered. Since the cross section falls down rapidly at 
increasing Q 2 , the absolute majority of the events occur at very low Q 2 , 
Q 2 = 0.5 GeV 2 . At such small Q 2 , the cross section a lN is theoretically 
proved to be a multiplier at inelastic fiA scattering sross section (Eq. (|6|)). 
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3 About possible physical reasons of rise of 
7 AT interaction cross section. 

Two problems arise on a way to reduce photon-hadron interactions to hadron- 
hadron ones. The first problem concerns the initial hadronic states. In 'jp 
interaction these states are a proton with the known quark-gluon structure 
and a quark bag, created by a photon. The structure of the latter is unknown 
and, moreover, statistical by its nature. In the framework of the AQM the 
problem is formulated as follows: what number of additive quarks describes 
the quark bag created by the photon? This question naturally arises, when 
one takes into account big widths of vector resonances. At transition of a pho- 
ton into p— meson, two alternative situations are possible: either p— meson 
is considered as an initial hadronic state and interacts with the proton as a 
single particle, or the initial state is the strongly bound tc + tt~ system having 



originated in the decay of p— meson [10]. One of the possible representations 



of this 7r + 7r~ system is a bag consisting of four additive quarks. In more 
general case the average number of additive quark-antiquark pairs must be 
considered function of energy of 7p interaction. 

The second problem is the dependence of photon hadronization proba- 
bility P 7 ^v r ( s ) on the energy. At ^Js > 10 GeV it can be higher than the 
standard value 1/250 on 10 — 20%, for example, due to the nondiagonal tran- 
sitions 7 — > V* — > V, where V* is a highly-excited resonance, V is a basic 
vector state (V = p, co, </)...) ]TU| . 



Let us turn to the phenomenological model that takes into account the 
above mentioned effects. Assume, that formation of the strongly bound 
7T + 7r _ system in the act of 7 — > V —>■ 2tt transition and its subsequent in- 
teraction with the proton is described by the model of four additive quarks. 
In this model the cross section of 2np interaction is estimated by the formula 



C27r P (s) = 2<Jn P (s/2) a np = -(cr-x+p + cr n - p ) 



1 

Let P v ->2ir{s) be the probability of conversion of a vector meson V into 2n- 
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system before interaction of V with the proton. Then the cross section of 7p 
interaction can be written as 

°7p( s ) = E P-y-*v{s) (1 - P V ->2n(s)) Ovp{s) + Pj->v(s)P V -^2ir(s)cr2Trp(s) 
v v 

= \ E Py^v( s ) (1 - ■ (3s/2) + 2 P 7 -,v(s)^-, 27 r(s)^(s/2) 



Formally substituting in (0) £P 7 ^ V = P^Xv = 1/250, P, 



(11) 

1, we 



obtain the so called conditional calibration curve 



4J(s) = 2P^ v W np (s/2) 



For its constructing we use data of groups JT3|, |14| . The fit of these data has 
resulted in the following expression 



(12) 



Ci = 1.155 /xb, > /iT= 1.672 GeV, A x = 152.0 /A, 

An = 5.006 GeV, A i2 = 6.658 GeV. 

Due to the universality of the Froissart asymptotic of hadronic cross sections, 
the parameter C\ = 2C is fixed by combinatorial reasons. In Fig.l the curve 
is represented by the upper solid line. 
Rewriting ([11]) in terms of calibration curves reads 



where 



a 7P (s) = k(s) \(l - Pv-ibr(s)) ^J(s) + P v ^ 2n (s)a$(s] 



Hs) = i + t^ot E AP 7 ^ y ^ y ( s ), 



p. 



(13) 
(14) 



AR^vis) = J2 P^v^v(s) - ifiv, Pi% * E P^v(s) 

V,V* V= P ,LO,<j> 

(see Eqs. ©, @). 
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The function k(s) takes into account the rise of photon hadronization 
probability due to nondiagonal transitions; P v ^2-k is the pionization proba- 
bility averaged over all vector mesons. 

One can easily see from Fig.l, that all measured cross sections of 7p 
interactions lie between the two calibration curves. This allows to assume, 
that the data can be interpreted within the framework of a hypothesis about 
the certain dependence of the pionization probability on energy or, what is 
the same, of the emergency probability quark bag consisting of four additive 
quarks. 

If we neglect all other factors increasing a photon hadronization proba- 
bility (in particular, the nondiagonal transitions) and put in (|13|) k(s) = 1, 
then the probability P v ^ 2 n{s) can be expressed only through the experimen- 
tal data: 

a - tr(°) 



rx (1) - rr (0) 



The values, extracted in such way and their fits by the formula 



P v ^2n(s) = wln n - (15) 

S2 



are shown in Fig. 2. 
The first fit, 



w = 0.0712, y/s2 = 5.465 GeV, n = 1/2 - for direct measurements; 

(16) 

represented by the dashed line, is based on the BUST data and the direct 
data [Rl ffl. The second fit, 



w 



0.0140, y/si = 2.166 GeV, n = 3/2 - for indirect data; (17) 



based on the BUST data and the indirect data 0, is shown by the dotted 
line. 

Substituting (IT5|) in fllBD, one derives the two curves for u w (s) fitting the 
two sets of experimental data mentioned above, which are presented in Fig. 
1. 
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Figure 2: Probability of vector meson pionization in 7p interaction. 

The results obtained, in our opinion, testify that the average number 
of additive quarks in the products of photon hadronization monotonously 
increases with energy under the law 



n^ h (s) = 2(l + w\n n ^j 



To obtain additional information, let us turn to the results of modeling of 77 
interaction cross section on the basis of the one. 
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4 Modeling of 77 interaction cross section and 
its comparison with the experimental data. 

The cross section of 77 interaction is expressed through the same probabilities 
of photon hadronization and cross sections of hadronic interactions, which 
appear in the model of 7p interactions in Section 3. The first step of the 
theoretical cross section calculation is obvious: 



<7 77 (sj 



1 — Pv-+2n-\ s , 



;is) 



Further on we express the cross sections of meson and two-pion interactions 
through the experimentally measured ones, simultaneously performing (ac- 
cording to the AQM) the scale transformation of arguments: 

a 77 ( S ) = k\s) (P(%) 2 (l - iW(s)) 2 ■ ^ <t Wp (9s/4) 

+2iW( S ) (l - • ~ a^s/A) + Pl^is) ■ \lf^/% ■ 

(19) 

At obtaining fll9"|) we use one more relation of the AQM: a wn (s) = cr^ p (3s/2)/cr Jvp (9s/4). 
The final result is convenient to present via the calibration curves: 



(7 77 y 3 



k\s) 



250 3 

+2P v ^ 2 Js) (l - P v -> 2 Js 



(l-P v ^As)) 2 -^(3s/2) 



a$(3s/2) + Pl^(s) ^ ' 



4?(3s/8) 



(20) 

In Fig. 3 the experimental data on er 77 (s) cross section from |27| and the 
theoretical predictions ( ^0|) at k(s) = 1 for two variants of photon hadroniza- 
tion probability flIT)| ) and ( |TTD are shown. 

The experimental uncertainties for this reaction, as well as those for 7^ 
one, do not allow to fix unambiguously the variant of high-energy behaviour 
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Figure 3: Cross section of 77 interaction and its modeling. 

of <7 7P (s) and cr 77 (s). However, the data on o" 77 (s) at y/s < 100 GeV, which 
have rather small errors, agree better with the forecast made on the basis of 
direct measurements of 7]? cross sections. 
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